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VIBRATION OF A WING OF FINITE SPAN IN A 
SUPERSONIC FLOW* 

By M. D. Haskind and S. V. Falkovich 

An investigation was made of the disturbed motion of a gas 
for the harmonic vibrations of a thin slightly cambered wing of 
finite span moving forward with supersonic velocity. This problem 
was considered by E. A. Krasilshchikova (reference l) who applied 
the method of Fourier series and obtained a solution of the space 
problem for the condition that the Mach cones drawn through the 
leading edge of the wing intersect the wing or are tangent to it. 

In this paper, a different method of solution is given, which 
is free from the previously mentioned condition. In particular, 
the vibrations of a triangular wing lying within the Mach cone are 
considered. 


1. FUNDAMENTAL EQUATIONS OF PROBLEM 

A thin, slightly cambered wing, whose projection on the 
x,y-plane has the form of an equilateral triangle with vertex 
angle 2y < 2a, where a is the Mach angle, is considered. It 
is assumed that the fundamental motion of the wing consists of a 
rectilinear translational motion with constant supersonic 
velocity u parallel to the x-axis . It Is also assumed that 
the coordinate axes move with the same velocity and that the 
x-axis is taken in the direction of the velocity u. On the 
fundamental motion of the wing, the additional harmonic of its 
vibration with frequency (0 is superposed, where the possi- 
bility of deformation of the wing is not excluded. The equation 
of the surface of the wing can then be written in the form 

z(x>y,t) s f Q (x,y) + f x (x,y) coscot + f 2 (x,y) sin cot (l.l) 


*"Kolebaniia Kryla Konechnogo Razmakha v Sverkhzvukovom Potoke." 
Prikladnaya Matematika I Mekhanika, Vol. XI, 1947, pp. 371-376. 
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It is assumed that the ratios of the functions to the lin- 

ear dimensions of the triangle and the derivatives bf-^/bx and 
bf k /bj are small . 

Considering the nonvortical motion of the gas, for the 
velocity potential $ (x,y,z,t) the following linearized equation 
may he taken; 

+ (i.2) 

3x^ 5y^ bz^ a xdt s? St^ 

where M « u/a is the Mach number and a is the velocity of sound. 

The flow condition that must he satisfied on the surface of 
the wing hy the velocity potential $>(x,y,z,t) shall he considered 
as satisfied on the projection of the wing in the plane z » 0. 

The derivation of this equation shall he based on the work of 
N. E. Kochin (reference 2). A stationary system of coordi- 
nates x-^, y-j_, and Zi connected with the moving system of 

coordinates hy the relations x = xi - uti, y » y^, z = zi, 
and t = ti is used j when expressed In these coordinates, equa- 
tion (1.1) "becomes 

Z 1 = f 0 ( x l- ut l^l) + 008 ^ + f 2 (*i-ttti,yi) sin COt 



For the component of the veloolty of the gas 
to the wing, there is obtained 

d z •» 3f q df -i Sf p 

— — m - U U r 003 <0t - COf n SlnCOt - U 

dt! Ox Ox bx 

By introducing the notation 

- u = Z(x,y) C0f 2 - u ■ Z 1 (x,y) -cof^ 

the boundary condition, which must he satisfied hy the velocity 
potential $(x,y,z,t), is represented in the form 


* Zo(x,y) + Zi(x,y) cos cot + Z2(x,y) sin cot (1.3) 
On the surface of the Mach cone, $(x,y,z,t) c 0. 



particles normal 


sin cot + cof^ cos cot 


bf. 


z 2 (*#y) 
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Thus the problem is arrived at: To find vithin the Mach cone 
the function $(x,y,z,t) satisfying equation (1.2) and condi- 
tion (1.3) and becoming zero on the Mach cone. 


When the steady vibrational character of the motion of the gas 
is taken into account, equation (1.2) is transformed to a simpler 
form, setting (reference 3) 


$> (x,y,z,t) = cp 0 (3:,y # z) + c Pi(**y> z ) cos (cat + Xx) + 



By introducing the notation 


Z 1 (x,y) + iZ 3 (x,y) «= Z(x,y) 
c P 1 (x,y,z) + icp 3 (x,y,z) = cp(x,y,z) 


there is obtained for determining the functions cp Q (x,y, z) and 
cp(x,y,z) after simple computations from equations (1.2) and (1.3) 


the following equations : 


d*tp A 
(l-M 2 ) -m 

ax 2 








ay 2 az 2 


= 0 


(: 


Sr)** - Z 0 (I > 7) 


d-M 2 ) + 1 * ♦ >*<?. o 

ax 2 ay 2 az 2 


> (1.4) 



i = iZ(z,y)e 


-iXx 


J 



(1.5) 


The function cp 0 (x,y,z) determined by equations (1.4) is the 
velocity potential of the motion of the gas for steady forward 

motion of the wing and Be [cp(3:,y,z)e’ i ^ t H is the velocity 

potential corresponding to the vibrations of the wing. 
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It is noted that for k = 0, the system (1.5) goes over into 
the system (1.4), so that it is sufficient to find the solution of 
system (1.5). 


2. METHOD OF SOLUTION 
In equations (1.5) 

X * X* VM 2 -1 

7 a 7* 
z = z* 


and the curvilinear system of coordinates are introduced 


\&2 " y * 2 -^* 2 
x*+y* 

z* 




x*+y* 


£ « ^/V* 2 -y* 2 -z* 2 


( 2 . 1 ) 




It is easily seen that within the Mach cone x* 2 -y* 2 -z* 2 = 0; 
the coordinates tj, and £ uniquely define the position of a point. 

On the Mach cone \ « 0, and £ = 0. In the triangular region of the 
projection of the wing, z * 0, |y|^ tg 7 x, 


7) = 0 

a ^ £ «T> 


A 




Voos 2y 


V'cos 27 


( 2 . 2 ) 


cos 7 + sin 7 cos 7 - sin 7^ 

In the new coordinates, equation (1.5) , assumes the form 


$2 ^ + 2 £ 2 / a 2 q> + sV 

at 2 at Vat 2 a n 2 / 


+ k 2 cp ■ 0 


(2.3) 


1264 
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'cW\ 

\Zt\/t\=0 


f(U) 


(2.4) 


where 


f(s,0 = Y iz (~i ] ~ 


£(1-£ 2 A / £(l+£ 2 )X r-p— \ 

— ; exp V ' 1 


By setting in equation (2.3) cp = X(| , T})\{/( £) after separation 
of the variables there results 


i 2 IF + 2t if + [* 2t2 - n(n+l % * - 0 

^ (if? + i^?) - n(n+1) x= 0 


(2.5) 

( 2 . 6 ) 


The solutions of equation (2.5) are expressed in Bessel functions. 
■When account is taken of the finiteness of the function cp on the 
Mach cone (£ = 0) , the following solution of equation (2.5) 
results : 


V'B ■ VS J «i (t£) < 2 - 7 > 

Y 2 

In order to choose in equations (2.5) and (2.6) the constant n, 
the arbitrary function F( £) is expanded, analytic in the region 
containing the origin of coordinates, into a series of 
functions 4r n ( £) : 

00 

F(£) - Y1 V'n^) 
n*0 


In order to establish the possibility of this expansion, use 
is made of the formula of Gegenbauer (reference 4, p. 283): 


i v _ 


W T >W*> (Itl < ItD 


( 2 . 8 ) 
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where A^t) ,are the Gegenbauer polynomials determined by the 
equation 


A nv (r) = L V ^.(V + n) ^ nu + n-m) n f » 

^o 1 m '- W 


(2.9) 


By setting u = | in equations (2.8) and (2.9), and by taking 
account of equation (2.7), the following series is obtained: 

OP 

fTf-I] V’WI) (l£!<|rp (2.10) 

* n=0 

where An(r) according to equation (2.9) will have the form 

'-n 


^2 n (1 \ 

vW*+l tr ? < ml \2 / 


( 2 . 11 ) 


VrtT n+1 m=0 

When series (2.10) and the integral formula of Cauchy are used. 


»(t) ■*&«> 

n=0 




«n 


2«1 


ff(7)A n (T)dT 


Tl <r 


(2.12) 


The coefficients can easily be expressed in terms of the 

derivatives of the function F( £) for £ » 0. Thus, by taking 
into account series (2.10) and the Cauchy formula for the 
derivatives, 

% , M^,( n -2i)( 0) (2n-2m)l 

2 n m^O ml (n-m) 1 (n-2m) 1 


(2.13) 


1264 
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For the steady motion in equation (2.5) k * 0, then 
\J/ n (£) = £ n and equations (2.12) hecome the usual power series 

no 

n*0 n* 

It is assumed that in equation (2.4) the function f(|,£) is 
analytic in then according to equations (2.12), it can he 

expanded in the functions ^ n (k£) into the series 


f(i A) - i f n ( W n (kt) (2.14) 

nsO 

In this connection, the solution of the problem is sought in 
the form 


<P« ,q,t) -y~!frnOct))Ue,n) (2.15) 

n=0 

On the basis of equation (2.4) and series (2.14), the boundary- 
condition for the function X n (l>T|), which satisfies equation (2.6), 
is obtained : 

(V 2 -) * f n (S) (a^b) (2.16) 

\Sri /ij-0 

where a and b are determined from equations (2.2). The problem 
thus reduces to the integration of equation (2.6) for the boundary- 
condition (2.16). By means of the transformation ^ *= | n+1 X n *, 
equation (2.6) reduces to the equation of D arbour for the 
function X n *(l,Tj) and because in this case n is an integer, 
the general solution of equation (2.6) can be represented in the 
form (see reference 5) 


^n( 


|n+l A n u n(£,n) 

t 


(2.17) 


where Uhdj^) I s 811 arbitrary harmonic function and A 11 is the 
Laplace operator of order n. 
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3 . OBTAINING THE HARMONIC FUNCTION 1% 

By developing the expression (2.17) and taking into account 
that AU n (! ,r|) * 0, there is obtained 


X n U,Tj) = (_2)m J2n^)i 


BteO 


ml (n-m) l 


gm 


c^Uh 


(3.1) 


When condition (2.16) is satisfied, an ordinary differential 
equation of Euler for (SUn/drOtyaO is obtained: 


f' , 2) m (|s^) i , B £L (?M . e a fnte ) 

n^O ml (n-m) l d| m \ Stj/tj=0 nl 


(as; |^b) 


(3.2) 


The general integral of an Euler equation, as is known, has the form 



(a=S Nb) (3.3) 


where !%(£) is a particular integral for equation (3.2) and C^' 
is an arbitrary constant. 


The complex variable T*= \ + itj is introduced and the 
function of this variable W n (T) * U n + iV n is considered. By 
making use of the Cauchy -Riemann equations, the boundary conditions 
that must be satisfied by the function W n (T ) are obtained: 


V n (|,0) = v n (|) = v n (£) °m£ 

DteO 




2m 




> (3.4) 


■U 


\jn(|) - u n(£) d £/ J 

Cm are new constants of integration and b . 


where 
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Farther, from equations (2.15) and (3.1), it follows that on 
the Mach cone 


<P(0,n,0) 


lim 






i n feo ' mi (n-m) l ’ d£ m 


When the properties of the function ^ n (k^) and equa- 
tions (2.1) are taken into account, 


cp(0 


»n»o) .£ ni(x+y) n UjO, 


n) 


n=0 


Because on the surface of the Mach cone, the function cp 
■becomes zero, it is necessary that U n (0,Tj) = 0. Moreover, on 
account of the symmetry 

Un(t,0) = 0 (0<| <a, |>b) (3.5) 

Condition (3.5) permits the analytical continuation of the 
function W n ( T) into the left half -plane. 

Finally, the function cp and therefore, according to equa- 
tion (3.1), the entire function ( t,q) must be finite for tj = 

and | j= a, and | = b. This condition will be satisfied if the 
finiteness of the function W n ( m )(T) is required for m = 0, 1, ...,i 
The satisfying of the last condition leads to equations determining 
the constants in equations (3.4). Thus the following boundary- 

condition problem of the theory of functions is arrived at; To 
determine the analytical function W n ( m )(T), regular in the upper 
half -plane, finite at the points T = ±a and t = ±b, and satis- 
fying the conditions on the real axis (tj = 0) : 

ImW n W - v n ( m )(£) ( +a< |<+b) 


Re W n ( m ) « 0 (a>[|| , j|j>b) 

ImW n ( m ) *= v n *(|) (-!><£< -a) 

(t»(-{) = (-1)" T n W({)) 
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By applying the formula of Keldish and Sedov (references 6 and 7), 


Wn (m, (T) 


W n W(T) 



STn^'fa) 


,2-t2 


v n (m) (a) 

s 2 -t 2 


/ b 2 - s 2 

a2_a.2 




ds (m-odd) 


>(3.6) 


s 2 -a 2 


ds (m-even) 


Equations (3.6) give for finite values at the points 

T * ia. The condition of finiteness of W n ' m )(r) at t =±b leads 
to the system of n+1 equations 



sM- v n ( m )(s) ds 


0 


(m*0,l, . . .,n) 


(3.7) 


where n = 0 for m even and p. « 1 for m odd. 

By computing v n ( m )(|) from equations (3.4) and substituting 
the result in equation (3.7), the system of equation (3.7) can be 
brought to the form 


Xt ° 


( Z®0, 1, . . . ,n) 


(3.8) 


where 


a 1 


2(Z-V) 


®Z,2V+1 “ 2(2-v)Bj^2v = 21(22-1) . . . (2Z-2V) , - Jg (l -p) 

i2U’V)dx 


* 


J 2(Z-1>) 


®Z 


nh 

s v v n (l) (s) ds 
a Y(s 2 -a 2 )(b 2 -s^) 


0 /\/(x 2 -l) (l-a 2 x 2 ) 

/ a f 0 for Z even \ 

^“b ' (ifor Z ' 


even 
odd J 
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t 

* Tli© integrals J 2 v can "be expressed in terms of complete 

elliptic integrals of* the first and second kind making use of the 
recurrence relation 


( 22+5 )k% ] /[_ ** 2(2+l)(l+k^)J2(2+iJ - (2Z+l)J2Z 


J 0 = Z* 


/v* \2 

J 2 « Z* + ( — ] E* 

\k/ 


(k« - tji-k 2 ) 


} (3.9) 


J 


where K* and E* are complete elliptic integrals of the first 
and second kind of the auxiliary modulus k* . All the coefficients 
Bj m of the system (3.8) 'can thus he computed. 


Translated hy S. Reiss* 
National Advisory Committee 
for Aeronautics . 
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